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Fig. 1. Convectively unstable sounding (with a moderately large midtropospheric 
minimum in the total static enthalpy), measured in the vicinity of a tornado at Jackson, 

MS on April 17, 1978 at 2300Z. The static temperature is T; the dew-point temperature, 

T<j; the pressure, p. 

Fig. 2. Pressure-temperature relations for the sounding of Fig. 1 (circles); for ground- 
level air raised on a moist adiabat to its level of neutral stability, at which its density is 
equal to that of the air of the sounding (triangles); and for neutral-stability-altitude air 
compressed dry-adiabatically to ground level (diamonds). 

Fig. 3. Schematic of a postulated four-part model of the structure of a quasisteady two- 
cell vortex, of axial extent zj (the “lid” of the storm). Region I is the outer potential vortex, 
II, the near-ground inflow layer; III, the “eyewall”, the lowest portion of which is termed 
the turnaround i and IV, the 4 *eye’\ The dotted curve is a sketched isobar, which is at altitude 

h at the periphery, r -> r e , and at lower altitude at the “eyewall’Vpotential-vortex interface, 
r = R(z), with R(0) = Ro and r e » Ro- The arrows schematically indicate the magnitude 
and direction of the secondary (i.e., radial and axial) flow. 


Fig 4 For the nominal case [defined by the sounding of Fig. 1 and the parameter 
assignments of Eq. (17)], the radial displacement of the eyewall, R(z), and the thickness 
of the eyewall, h(z). 

Fig. 5. For the nominal case, the streamwise velocity component in the eyewall, q(z), 
and the swirl velocity component, v(z). 

Fig. 6. For the nominal case, the pressure at the potential-vortex/eyewall interface, 
Pamb( z )' and the height at the periphery at which the pressure pamb occurs, h(z). 

Fig. 7. For the nominal case, the mass in the eye as a function of altitude, meye(z), 
and the filling time T(z)= m e ye(z)/rh e yewaii, where rh e y e waii is defined by Eq. (14). 


Fig. 8. The radial displacement of the eyewall, R(z), for several values of the parameter p, 
related to the vertical distribution of the angular momentum at the periphery, where p = 1 
is an axially invariant distribution and p = 0 is strongly stratified. All other parameters and 
functions are held at their nominal values (and p = 0.99 is nominal). 42 

Fig. 9. The eyewall thickness, h(z), to complement the results of Fig. 8. 43 

Fig. 10. The radial displacement of the eyewall, R(z), for three values of the initial thickness ho, 
all other parameters and functions being held at their nominal values (with ho — . 

nominal). Of course, the initial curvature e [Eq. (18)] varies with ho, since the ratio (Rc/ho) is 
fixed at its nominal value, 10. 


Fig. 1 1. The eyewall thickness, h(z), to complement the results of Fig. 10. 
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Fig. 12. A simplified schematic of the structure of an axisymmetric two-cell vortex 
convenient for thermohydrostatic computation. The surface-friction layer II and 
the turnaround (Fig. 3) have been omitted; in the above sketch, the surface frictional layer 
does not separate smoothly, and the pressure at the base of the eyewall (region HI) 
appreciably exceeds that at the base of the eye (region IV). For a less intense, one-cell 
vortex, no eye is present, and r 0 -» 0, so region HI (no longer termed an eyewall) 
envelopes the axis of symmetry. 

Fig. 13. The total-static-enthalpy profile H(z) associated with three columns of fluid for 
the model of Fig. 12, with the sounding of Fig. 1 holding at the periphery. 

Fig. 14. The pressure profile p(z) associated with three columns of fluid for the model of 
Fig. 12, with the sounding of Fig. 1 holding at the periphery. 


Fig. 15. The temperature profile T(z) associated with three columns of fluid for the 
model of Fig. 12, with the sounding of Fig. 1 holding at the periphery. 

Fig. 16. The density profile p(z) associated with three columns of fluid for the model of 
Fig. 12, with the sounding of Fig. 1 holding at the periphery. 


Fig. 17. For spatially constant viscosity, results related to the radial and axial velocity 
components for the outer, preponderant portion of the boundary layer under the high-speed 
portion of an impressed vortex. The similarly independent variable r\ is large for large 
distances normal to the wall or for small radial distances from the outer edge of the 
vortex system. 

Fig. 18. As a complement to Fig. 17, results related to the swirl-velocity component. 

Fig. 19. For spatially constant viscosity, a comparison of numerically computed (dotted 
curve) and approximate analytic (solid) results for the function h(a), related to the radial 
inflow in the thin, effectively nonrotating sublayer under the high-speed portion of an 
impressed swirl. The similarity independent variable o is large for large distances normal 
to the wall or for small radial distances from the axis of rotation. 


Fig. 20. For numerical integration of the turnaround, a convenient (but nonessential) eight- 
segment subdivision is considered, with a segment end-point introduced wherever 

zr = 0, Rz = 0, or (Rzl= 1. Equation I is Eq. (D.8); Eq. II is Eq. (D.12), derivable from Eq. (D.8) 
by an interchange of the roles of the independent and dependent variables. The “sign” 
column indicates the sign locally adopted for the last turn in Eq. (D.8) or Eq. (D.12), as 
appropriate. ^ 

Fig. 21. The dimensionless radial displacement, R(z), of a smoothly separated surface- 
inflow layer, over one cycle (incursion, then excursion, with increasingattitude z) of the 
turnaround. As the geometric ratio a increases, the minimum value of R increases, so the 
local peak swirl exceeds the thermohydrostatically estimated peak value by a smaller 
amount. ^5 
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1. INTRODUCTION 


We adopt an idealized model of an intense tall axisymmetric vortex in a stratified “moist 
atmosphere” above a ground plane, and then analyze the flow to ascertain whether the vortex could 
be selfsustaining. The ultimate objective is the identification of (hopefully satellite-accessible) 
observables for anticipating tomadogenesis within a supercell thunderstorm. 

In undertaking this work, we note the assessment of Rotunno (1986) that, while aircraft and 
multi-Doppler radar with a resolution of few kilometers have increased understanding on the 
10-50-km thunderstorm scale, the more difficult question of how the larger tornado cyclone 
produces the smaller- scale tornado is unresolved. “Given the present state of computers, it is 
impossible for three-dimensional cloud models, even with nested grids, to simulate both the 
supercell on a domain large enough to include the undisturbed environment and the tornado with a 
resolution fine enough to capture the significant dynamical processes. However, even if this were 
possible, the results of such a calculation might be too difficult to handle.” (Rotunno 1986, 
p. 427.) 

Supercell-spawned tornadoes have the highest wind speeds, the longest and widest paths, and 
the longest life-span of any tornadoes, and, thus, from a practical point of view, are the most 
important (Rotunno 1986). We know that severe weather arises in convectively unstable ambients 
(i.e., tropospheric stratifications in which the low-level air is particularly warm and moist, so there 
is a pronounced midtropospheric minimum in the total static enthalpy) (Darkow 1986). However, 
convective instability only sometimes results in a supercell thunderstorm, and perhaps only one in 
four supercell thunderstorms spawns a tornado. Thus, a tornado warnings based on mesocyclone 
detection by dual-Doppler radar would seem likely to incur a significant false-alarm rate, although 
mesocyclones often manifest some type of severe weather (Brown 1987). 

Identifying a further discriminant for tomadogenesis in already well-developed, rotating, 
severe local storms remains a worthwhile goal. Inquiring into the vertical-vorticity-amplification 
mechanisms responsible for the pretomadic organization of the rotating thunderstorm seems 
nonessential for this objective; the organization may be viewed as a given “initial condition”. 

Rotunno notes that, in the present state of limited understanding of the asymmetric 
heterogeneous flow in the storm, guidance often is sought from examining the axisymmetric 
homogeneous flow in a cylindrical^ configured laboratory apparatus called a vortex chamber 
(Maxworthy 1982). Fluid within such a chamber is drawn through it by suction applied at an 
exhaust hole in one end wall, and rotation is added to new fluid, which enters through the lateral 
surface of the chamber. Clearly, both the primary flow (circulation) and the secondary flow 
(inflow-updraft-outflow) are very highly constrained in a vortex chamber, and the circulation 
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added to entering air typically does not vary with height Thermodynamically, there is no 
condensation of water vapor, and no density stratification with altitude because of the limited 
height of the chamber. The laboratory vortex chamber may be of value for investigating specific, 
purely fluid-dynamical aspects, such as the mass exchange between contiguous “columns” of air 
(with differential rotation) in a vortex, although artifacts of the apparatus could be misleading about 
other properties of an intense atmospheric vortex. In our concluding remarks, we suggest a further 
aspect of tornado dynamics that may be usefully pursued in a laboratory vortex chamber. 

For elucidating the thermohydrodynamics of tomadogenesis, we adopt a model of rapidly 
swirling flow that: extends over the vertical scale of the troposphere; is axisymmetric (for analytic 
convenience); and involves moist air (with both static temperature and dew-point temperature 
stratified in accord with atmospheric soundings) which undergoes idealized processes. In selecting 
appropriate idealized processes, we note that the inferred peak swirl speed, relative to the axis of 
rotation, in the exceptional, intense tornado (of interest here) may exceed 100 m/s. We believe 
that the only plausible tropospheric mechanism for attaining the magnitude of lateral pressure 
deficit (near ground level) consistent with such high swirl speed is the “two-cell”, “eye”-within-an- 
“eyewall” structure known to be associated with hurricane- stage tropical cyclones (Fendell 1974). 
We suggest that some phenomena associated with the intensification of a tropical storm into a 
hurricane-stage tropical cyclone also occur during tomadogenesis in a mesocyclone, albeit on 
vastly reduced lateral scale (Dergarabedian and Fendell 1970, 1971a, 1971b, 1977, Walko 1988). 
Specifically, we associate tomadogenesis with the (at least partial) insertion of a central column of 
dry, only very slowly rotating, very slowly recirculating air, characterized by a nearly dry-adiabatic 
locus of thermodynamic states, within an annulus of mostly saturated, rapidly rotating and 
ascending air, characterized by a nearly moist-adiabatic locus of thermodynamic states (Figs. 1 and 
2). Of course, the idealization of a nearly dry “eye” and an entrainment-free “eyewall” is a 
convenient simplification that assists both the exposition and the analysis below, relaxation from 
this idealization may be the subject of subsequent studies. The compatibility of exceptional swirl 
speeds with this just-described two-cell structure has been previously reported by considering 
measured, convectively unstable ambients under simplifying cyclostrophic and hydrostatic 
approximations. 


2. OBJECTIVE AND MODEL 

Our objective is to examine a further point about the above-described swirling flow by 
applying the laws of dynamics and thermodynamics. We inquire whether only certain 
stratifications of the ambient angular momentum at the periphery of the flow permit the two-cell 
structure, requisite for intense swirl speeds of roughly 100 m/s, to persist. For example, perhaps 
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a two-cell structure is plausibly sustainable only if the low-level angular momentum at the 
perimeter persists on a scale equal to the height of the troposphere. 

We adopt the model in Fig. 3 for investigating the conjecture of the last paragraph. Although 
we undertake no transient analysis, we believe that the adopted four-part quasisteady structure 
might well evolve after spin-up from “reasonable” initial conditions. The top of the flow (say, near 
15 km or so) is defined to be an isothermal isobaric horizontal plane, located at an altitude at which 
the value of the ground-level total static enthalpy is virtually recovered, for the convectively 
unstable stratification holding at the “lateral edge” (or periphery) of the flow field. At the 
periphery, the static temperature, the dew-point temperature, and the angular momentum are 
explicitly given (or are readily deducible) as functions of altitude above the ground plane. (As 
noted, at the periphery, we typically shall adopt temperature and dew-point-temperature profiles 
measured in temporal and spatial proximity to a tomadic thunderstorm.) The ground plane is an 
impervious isothermal plane at which the no-slip boundary condition holds. In the next paragraph, 
we provide more detail on the quasisteady four-part structure, depicted in Fig. 3, much of the 
remainder of this investigation then seeks to identify the ambient conditions under which the four- 
part structure is plausible, stable, and compatible with thermohydrodynamic principles. 

Most of the flow (except near the ground plane and the axis of rotation) would consist of 
warm moist air rotating in a potential vortex in a cyclostrophic, hydrostatic balance, such that 
isobars would dip groundward as the (cylindrical) radial distance from the axis of symmetry (and 
rotation) decreases (region I of Fig. 3). There is very little radial motion in the spun-up potential 
vortex, and only a relatively slow downflux into a nonlinear frictional layer (region II), contiguous 
to the ground plane. The frictional-layer thickness increases modestly (to, say, 50 m or so) with 
decreasing radial distance from the axis. Rotating flux that descends into this frictional layer flows 
radially inward within the layer; in fact, for an intense vortex, al smaller radial distances from the 
axis of symmetry, frictional effects are significant only in an ever-thinning sublayer immediately 
contiguous to the ground plane, and most of the inflow-layer thickness can be described by 
inviscid dynamics. The swirling inflow layer separates at a finite radial distance (say, at radius of 
50-75 m) to become a swirling annular updraft of saturated air (the “eyewall”, region HI). The 
very-low-altitude comer-like flow of region in is henceforth referred to as the turnaround, and is 
treated separately below by specializing the relations holding more generally over the vertical extent 
of the eyewall. In Fig. 3, the eyewall becomes thicker, and tilts away from the axis, with 
increasing altitude; the updraft may well be sufficiently vigorous that little condensate falls out 
Within the updraft annulus is a central column (the “eye”, region IV) of warm, dry, slowly 
recirculating, only very slowly rotating air. For a finite amount of warm moist air available for 
throughput, the swirling effectively dies (or at least is very diminished) after all the outer-vortex 
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fluid has been “processed”, owing to the frictional loss of angular momentum in the near-ground- 
level-inflow layer. 


3. A MODEL FOR THE ANNULAR ROTATING UPDRAFT 

A transversely averaged model for the “eyewall” (region HI) is given by the following set of 
coupled algebraic and differential equations [previously presented in more cumbersome form and 
without solution (Carrier et al. 1984)]: 

conservation of mass (with entrainment from the eye into the eyewall), 

p q R h = p 0 q 0 Ro ho + a E I p ey e(zi) q(zi)[R(zi) * Kzi)] dzi ; 

Jo ( 1 ) 


conservation of angular momentum, 


Rv = Ro v 0 ; 

conservation of radial momentum across the annulus, 

P<l 2h d 2 R P v2 h 


[l + (dR/dz) 2 ] 3/2 dz2 R[l -h ( dR/dzp] 1/2 


+ Pamb “ Peye — 0 » 


( 2 ) 

(3) 


Bernoulli's equation. 


q 2 + v 2 

2 


+ gz + 



dp q 0 2 + v» 2 . 

P(P) 2 ’ 


(4) 


moist-adiabatic locus of thermodynamic states based on the ground-level state holding at the 
periphery (previously tabulated and presented here formally), 

p = F(p) ; (5) 

definition of a "mean eyewall pressure", 

2p = Pamb + Peye * (6) 

and selfconsistency of an isobar in the outer, potential vortex [see Eqs. (1 1)-(13) below and 
Appendix C], 

h(p,mb) = [ t i ■ (1 ■ bz)] exp|- b j . (7) 
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These seven equations are to be solved to yield the dependent variables q, v, p, p, Pamb> R. and h 
as functions of the independent spatial variable z(> 0), with the starting conditions at z = 0: 

R(0) = R 0 , = N , given , (8) 

where we expect N = 0 to be of particular interest. That is, the above equations are taken to hold 
for an initially vertical updraft, after the throughput has transversed the turnaround, and possibly 
undergone a “vortex breakdown”, in the very-low-altitude subdomain of region IE. As yet, we 
have not studied the selfconsistency of the conditions of Eq. (8) with the behavior holding in the 
turnaround. In any case, the origin z = 0 is still “close” to ground level. Subscript o denotes 
known values holding at z = 0. The function p ey e(z). which denotes the pressure as a function of 
altitude z in the central dry, transversely invariant core, is known from preliminary 
“thermohydrostatic” calculation, as is the tabulated, moist-adiabat-based relationship [Eq. (5)] 
between p and p, the density and pressure, respectively, in the updraft annulus. For 
completeness, this preliminary set of calculations is summarized in Appendix A. 

The (cylindrical-)radial displacement of the “eyewall” from the axis of rotation (and symmetry) 
is denoted R; the thickness of the “eyewall”, h; the swirl speed in the “eyewall , v; the flow speed 
along the “eyewall”, q; the pressure at the interface between the “eyewall” and the outer potential 
vortex, pamb; and the magnitude of the gravitational acceleration, g. Further, from Appendix A, 

V 0 = [2{Ap)/po] m . (9) 

where (Ap) is the near-ground lateral pressure deficit (between the axis of symmetry and the 
periphery) for a two-cell structure. Also, if q G is the average radial velocity in the boundary-layer 
flux, and v 0 is the average swirl velocity in the boundary-layer flux, 

V? = qg + vi , (10) 

a relation obtained in the course of the analysis of the higher-peak-speed, nearer-to-the-axis-of- 
symmetry portion of the surface frictional layer (Appendix B), in which it is suggested that q 0 - v 0 
in the frictional layer near separation. However, q 0 (but not v G ) may be decremented in any 
subsequent vortex breakdown, and the decremented value of q G would be appropriate for use in 

Eqs. (1) and (4). 

As developed in Appendix C, the swirl in the outer potential vortex is given by [r is recalled to 
be the (cylindrical-)radial coordinate] 

v|u,„= Y ^r sif ( z ). an 
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where at the outset we study the specific case of a linear distribution of the peripheral angular 
momentum with altitude in the troposphere: 

Hz)=P + (1-P)(^).0SP<1. < 12 > 

We let 

f(z) s a - bz , (13) 

so a = 1, b = (1 - P)/zt- That is, (5 = 0 implies a peripheral swirl so stratified with altitude that no 
swirl exists at the "lid" z = zt, the known height of the "tropopause"; alternatively, p = 1 implies a 
distribution of swirl invariant with altitude. While variation of results with the parameter P is of 
key interest, we note here, and reiterate below, that if the swirl in the eyewall is smaller that the 
swirl in the potential vortex near the potential-vortex/eyewall interface, the flow is unstable. 

Hence, a selfsustaining vortex would not be expected for smaller values of the parameter p. 

Eq. (1 1) links the altitude h, at which a particular value of the pressure Pamb arises at the 
periphery, with the (lower) height z at which that same particular value of the pressure pamb holds 
at the "eyewalTVpotential-vortex interface. The function h(p), which gives the height associated 
with pressure at the periphery, for the peripheral sounding under investigation, is obtained, of 
course, by hydrostatics from the given data (static temperature and dew-point temperature as a 
function of pressure). 

A factor of 2n has been cancelled from each side of Eq. (1), in which the Taylor entrainment 
constant oce - 0.08 - 0.09 (Turner 1969). While we anticipate entrainment of dry, warm, virtually 
nonrotating air from the eye into the eyewall, we do not expect much entrainment into the eyewall 
from the contiguous portions of the potential vortex, especially for those altitudes for which the 
contiguous-potential-vortex fluid is rotating at a speed in excess of the transversely averaged swirl 
speed holding in the eyewall. Since the mass entrained into the eyewall is taken to be detrained 
from the eye, an axial downdrift w eye arises in the eye. If the eye/eyewall interface at altitude z 
occurs at (cylindrical-)radial distance [R(z) - h(z)], we may estimate the downdrift by the 
conservation of mass: 

p eye (z) [R(z) - h(z)] 2 W eye (z) = -2(Xe j Peyefcl) q( z l)[R( z l) ' K z l)l <* z l » 

where we have taken the downdrift to vanish at the impervious ground plane z = 0. We anticipate 
the downdrift to be comparatively modest, for consistency with the slow recirculatory flow 
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envisioned for the eye relative to the rapid updraft in the eyewall. The influx into the eye in the 
upper troposphere is not within the scope of this discussion. 

The appearance of the two principal radii of curvature in Eq. (3) may be noted. Equation (3) 
is to be integrated by: (1) holding the product (ph) and the difference (pamb - Peye) constant over a 
very small increment in altitude z; (2) using the modified value thus obtained for the eyewall 
radius, R, to update the values for p, h, p eye > and pamb at the slightly incremented altitude, and 
(3) repeating steps (1) and (2) for a further, small increment in altitude, until one attains either the 
highest altitude, zj, for which the thermodynamic and dynamic state at the periphery is defined 
(Appendix A), or until singular behavior intrudes. 

We anticipate that, in general, the slope of the eyewall, R z , is small. If not, then it becomes 
noteworthy that, whereas all the other dependent variables are functions of the altitude z (as 
previously stated), in fact the argument of the tabulated function p eye is z', where 


For completeness, we note that we have regarded the eyewall/potential-vortex interface to be 
the definition of the displacement R(z), for circumstances in which the nominal smallness of the 
ratio h/R does not hold. 

While the spatial scales Ro and h G may appear to be independently assignable, in fact, a thin- 
layer-type, transversely averaged treatment of the eyewall is more appropriate for circumstances for 
which the ratio (Ro/h 0 ) > 0(10). More specifically, the case (R</ho) = °( 10 ) is more P lausible 
under the observation that the mass flux entering the eyewall, m, is given by 

rileyewall 5=5 2ft Po ^0 Qo * (14) 


the mass content of the eye, m e ye> is given by 


m e ye — 


J .*r f R u) 

*1 1 

0 Jo 


Peye( z ) 2ft T 


„..r 

Jo 


R^(z) p e ye(z) dz , 


(15) 


where p eye (z) is available from preliminary thermohydrostatic calculation (Appendix A), and R(z) 
is obtained from solution of the problem posed for the eyewall. If we envision the filling of the 
eye from throughput detrained from the eyewall during a spin-up stage, which occurs pnor to the 
quasisteady mature stage under examination, then the time T required for such filling is given by 

T — m e ye /m eyewall • (10 
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The upper limit zj in the integral for nieye may be taken to be a somewhat lower value if, as z 
approaches zj, the eyewall displacement R slopes appreciably away from the axis of symmetry. 
(The eye/eyewall structure probably needs re-examination in the form of an outflow-layer analysis 
for the portion of the flow near the lid.) However, we anticipate that we shall find that the time for 
filling, T, exceeds the observed lifespan of the tornado if Rq far exceeds ho; if so, such a 
configuration seems implausible. 
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4. RESULTS FOR THE ANNULAR ROTATING UPDRAFT 


With a preliminary treatment in hand for the eyewall formulation (Section 3) specialized for the 
turnaround (Appendix D), we now seek results to Eqs. (1)-(16) for the bulk of region III (Fig. 3). 

For brevity of reference, we henceforth define the following set of parametric and functional 
assignments as “the nominal case.” The thermohydrostatically derived input z T > V 0 , Peye(z). Po. 
and p(p)-the density-pressure relation for the eyewall-is based on the sounding given in Fig. 1 ; 
for convenience, we note that z T - 1 1.6 km, V D - 106.9 m/s, Po - L12 kg/m*, Peye(z) is given in 
Fig. 14, and p(p) is given parametrically in terms of the altitude z by Figs. 14 and 16. Since the 
gravitational acceleration g = 9.8 m/s 2 and v D = V 0 /2 1/2 , the remaining quantities to be assigned 
are p, a dimensionless parameter related to vertical stratification of the peripheral angular 
momentum [see Eqs. (1 1)-(13)]; M, a dimensionless parameter related to kinetic-energy loss in the 
breakdown [see Eq. (B.36)]; N, a dimensionless parameter related to initial eyewall slope [see 
Eq. (8)]; h G , the starting thickness of the annulus; and Ro, the starting radial displacement of the 
eyewall (or equivalent data). For the nominal case, we choose 

ot E = o, P =s 0.99, M = 2 1/2 , N = 0, ho = 35 m, and (Ro/ho) * 10 , (17) 

where we discuss implications of the last assignment in the next paragraph. Operationally, we find 
that an integration step size Az - 0.5 m yields results that are effectively invariant with further 
refinement of the discretization; with Az - 0.5 m, solving Eqs. (1)-(16) over 0 < z < z T requires 
roughly 15 min on a 486 PC. 

If we let e = R zz (0) for brevity of notation, then at z = 0, from Eq. (3), 


£ho = 



PqVq/Sp 

Ro/ho. 



(18a) 


PoQ^/Sp 


p 0 v o/5p _ 1 


Ro / ho 


(18b) 


for the case of primary interest, N = 0 ; 8p s Pamb (0) - p ey e(0). Thus, in general, a finite positive 
(but typically small) curvature exists at z = 0. By rearrangement of Eq. (18b), 


R„ Po v o /S P 

ho ehdjp 0 q£/8p) + 1 (l gc ) 
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We see that a smaller (but still positive) value of the nondimensionalized curvature, eho, implies a 
larger value of the ratio (Rq / ho); in fact, eho = 0 typically implies an appreciably larger value of 
(Ro/ho) than that given in Eq. (17), with possible further implications of implausible tomadic 
structure [recall Eqs. (14) and (15)]. 

For the nominal case. Fig. 4 presents the radial displacement R(z) and the annular thickness 
h(z); Fig. 5, the streamwise velocity component q and the swirl v; Fig. 6, the pressure at the 
potential-vortex/eyewall interface, Pamb( z )* and the height at the periphery, h(z), at which the 
pressure pamb( z ) occurs; and Fig. 7, the eye mass m eye (z) and the filling time T, where the upper 
limit of the integral Eq. (15) has been changed from zj to the general value z. The value of the 
difference [pamb( z ) * Peye( z )l is quite small relative to the value of either contribution, and varies 
from about 500 Pa near an altitude of 1 km to about 300 Pa at 3 km to about 100 Pa at 6 km; the 
difference goes to zero at about 8.5 km and even very slightly negative above that altitude. The 
filling time at z = 8.5 km is about 20 min, but increases rapidly for higher z because the 
displacement R(z) becomes large at large values of z. The annular thickness h monotonically 
increases from 35 m at z = 0 to about 64 m at 5.8 km, and then monotonically decreases for greater 
height; correspondingly, the streamwise speed q decreases from 75.6 m/s at z = 0 to about 38 m/s 
at 6.8 km, and then monotonically increases for greater height 

Figures 8 and 10 presents R(z) and h(z), respectively, for several values of the parameter p, all 
other parameters and functions being held at their nominal values. For the more strongly stratified 
profiles of the peripheral angular momentum, the eyewall remains very vertical. While a solution 
may be obtained formally for the strongly stratified case p -* 0, its stability warrants investigation. 
Stability requires that the angular momentum should not decrease with radius, at fixed attitude. 
Thus the solution is unstable if the inequality v(z) < v outer [R(z), z] holds for any z, where v oute r is 
defined by Eq. (11). 

Figures 10 and 1 1 present R(z) and h(z), respectively, for three values of the parameter ho, 
with all other parameters and functions fixed at their nominal values [e.g., (Ro/ho) = 10]. Thus, 
for the nominal case h 0 = 35 m, Rq « 350 m (and £ = 0.002); for h Q = 17.5 m, Rq *= 175 m (and 
£ « 0.004); and for h Q = 70 m, Ro » 700 m (and £ = 0.001). Accordingly, the mass flux 
p 0 h 0 Ro q Q varies by a factor of 16, and the product Ro v 0 by a factor of 4, over the three cases. 
The time T to fill the eye effectively does not vary for these three cases. 

Further parametric investigation with the analytic framework already provided, and with the 
extensions briefly noted in the next section, should permit additional insight into criteria for the 
existence of selfsustaining intense vortices. 
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5. CONCLUDING REMARKS 


The observable under scrutiny as a possible predictor of tomadogenesis in a raesocyclone is 
the vertical profile of the angular momentum at the periphery of the mesocyclone. This quantity is 
closely related to the vertical profile of the circulation, and if indeed the rotation within a tomadic 
mesocyclone is very roughly that of a potential vortex (except in the core), then the vertical profile 
of the circulation should be invariant with the lateral distance from the axis at which the profile is 
measured, for sufficiently large lateral distance. A single sounding is inadequate to obtain the 
circulation as a function of altitude in the vicinity of a mesocyclone in the atmosphere, and 
probably futuristic, continuously monitoring, geosynchronous-meteorological-satellite-based wind 

sensors with vertical resolution would ultimately be required. 

Among the future improvements and generalizations of the model is the desirability of 
smoothing the transition of the total static enthalpy H at the interface of the eyewall (within which 
H(z) = Heo) and the potential vortex [with which H(z) = He(z)]. 

Still another topic warranting further investigation is the upper-tropospheric radial outflow of 
“processed throughput.” In a closed system, this throughput would be compressionally heated as 
it descended to accommodate the further processing of the convectively unstable fluid of the 
potential vortex; while evaporative cooling afforded by small-sized condensate carried along by the 
flow would somewhat counteract the compressional heating, still the descent would partly negate 
the lightening effected by the presence of an eye. In fact, we believe that there is a high-level 
lateral efflux of processed air. 

However, we particularly emphasize the desirability of seeking further insight into turnaround 
dynamics by use of a laboratory vortex chamber. A theoretical investigation of the complicated 
turbulent flow of a compressible fluid in the unconventional geometry of the turnaround would be 
assisted by experimental verification. Thermodynamics is not at issue in this particular region; 
inability at laboratory scales to experimentally imitate the moist-air behavior of an atmospheric 
storm is irrelevant 
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APPENDIX A. THERMOHYDROSTATICS FOR A ONE-CELL OR TWO-CELL VORTEX IN 
CONVECnVELY-UNSTABLY -STRATIFIED MOIST AIR 

A.l Objectives and Idealizations 

By use of (1) the thermodynamics of moist air, (2) the conservation of vertical momentum in 
the hydrostatic approximation, (3) a sounding (of the static temperature and the dew-point 
temperature as a function of pressure for a convectively unstable atmosphere), and (4) a simplistic 
axisymmetric quasisteady model of a one-cell vortex and of a two-cell vortex (Fig. 12), we obtain 
the thermodynamic states holding within the flow. For a one-cell vortex, fluid at and near the axis 
of symmetry is described by the locus of thermodynamic states associated with a moist adiabat, 
based on the ground-level state at the peripheiy. The associated low-altitude lateral pressure deficit 
from axis to periphery is relatively modest. For a two-cell vortex, the fluid ascending (nearly 
vertically) on a moist-adiabatic locus rises in an annulus at a finite radial distance from the axis; the 
fluid in region IV of Fig. 3 (or Fig. 12) has a locus of thermodynamic states associated with a dry 
adiabat, based on the thermodynamic state holding at the top of the flow (at all radial distance from 
axis to periphery). The eye-insertion-associated low-altitude pressure deficit from axis to 
periphery is relatively large. Since perfectly entrainment-free ascent on a moist adiabat is an 
idealization, as is slowly recirculatory flow in a perfectly dry eye extending from the top to the 
bottom of the vortex, the pressure deficits computed for such idealized processes are upper bounds 
on the corresponding pressure deficits that would be observed for real counterparts of the 
processes. 

The low-altitude peak swirl speed associated with the just-discussed axis-to-periphery 
pressure deficits may be computed by use of (1) the cyclostrophic approximation to the 
conservation of angular momentum; (2) an adopted dependence of the swirl on the radial 
coordinate [the dependence being taken to be a Rankine (combined) vortex for one-cell structure, 
and a rotationless core patched to a potential vortex for two-cell structure]; and (3) a radial profile 
for the density [taken to be constant at the ground-level value for the moist- adiabatic column]. As 
just noted, the lateral pressure deficit is appreciably smaller for a one-cell vortex, and half of the 
pressure deficit is expended to maintain the rigidly rotating core of a Rankine vortex; thus, the 
thermohydrostatically derived peak swirl speed for a one-cell vortex is relatively modest, say, 

50 m/s or less. The lateral pressure deficit is appreciably greater for a two-cell vortex, and no 
portion of the deficit need be expended to support rotation of the core, which is virtually 
nonrotating; thus, the thermohydrostatically derived peak swirl speed for a two-cell vortex is 
relatively large, say, 100 m/s or more. 
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There is an altitude at which that total static enthalpy of the near-surface air at the periphery is 
recovered in the upper troposphere. In fact, we may refer to the height of the lid as the tropopause, 
though this is by no means the universally standard definition. 

A.2 Collection of Standard Relations; Introduction of Notation 

If subscripts a and v refer to dry air and water vapor, respectively, then, for convenience of 
future reference, we recall that the pressure p and the density p obey 

p = Pa + p v ,p = Pa + Pv- (A.1) 

Also, 

Y a = ^,Y v s^,soY a + Y v =l (A.2) 

in the absence of condensed water substance. The equations of state for dry air and water vapor, 
respectively, if R denotes the gas constant for dry air and a denotes the ratio of the molecular 
weight of water vapor to that of dry air, 

p a = p a RT, ap v = pv RT, so p - (1 - a) p v = p RT , (A.3a) 


and 

Y v = cp v /[p -(!-<*) Pv] • (A. 3b) 

The saturation vapor pressure for water is denoted P(T), where a standard excellent curve fit 
to data (Murray 1967) is given by, with P in Pascals and T in Kelvin, 


P(T) = 6.1078 x 10 2 exp[a(T - 273.16)/(T - b)] , (A.4a) 


where 


a 

b 


21.8745585 

7.66 


over ice. 


a = 17.2693882 
b = 35.86 


over water . 


(A.4b) 
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By definition, the relative humidity RH and the dew-point temperature T<j are related by the 
following: 


RH = P y = 

P(T) P(T) * 

Thus, if the static temperature is given, RH is implied by T<j, and vice versa. 
The total static enthalpy H, by definition, is 


H = CpT + gz + LY v , 

where, in standard notation, 


(A.5) 


(A.6) 


c p - Y a Cp, + Y v Cpv ~ Cp, , (A.7) 

since (Y v /Ya) « 1 and c pi /c pv = 0(1) for all the circumstances of interest. [Incidentally, at ground 
level in warm moist air, T = 300 K and (L Y v /c p ) = 40 K; for comparison, if we introduced a total 
enthalpy (including dynamics), the typical peak tomadic-kinetic-energy contribution 
Vo / (2 c p ) = 4 K.]. We also ignore the variation of the specific latent heat L with temperature, and 
take all condensed water vapor to remain as liquid water (so no heat of freezing is realized). 

Finally we recall that, under hydrostatics, 

3p 

— = - Pg • (A.8) 


so Eq. (A.6) may be written as 


dH = c p dT - ^ + L dY v 
P 


(A.9) 


For a constant-static-enthalpy process in unsaturated air, so virtually no condensation occurs and 
dY v = 0, then, with Y v = (Y v )h;f, const.. 


dT _ 1 _RTf 
dp pc p P l 


1 + 



(A. 10a) 


or, with (R7c p ) = (y - l)/y, 


T 

Tref 



(A. 10b) 
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The other thermodynamic variables p, p v , pa» etc., follow directly. Alternatively, for a constant- 
static-enthalpy process in saturated air, so RH « 1, and p v = P(T), 

dY v = (c[p'(T)] p dT - a P(T) dp) / x x s p - (1 - a) P(T) , (A.l 1) 

and 


RT , LgP(T) 
dT _ * x2 

dp lopttip' 

t : 


(A. 12) 


Again, the other thermodynamic variables follow by straightforward substitution. For (Y v )tef = 0, 
Eq. (A. 10b) describes the locus of thermodynamic states on a “dry adiabat”. The locus of states 
on a “moist adiabat” is described by Eq. (A. 10b) prior to saturation (RH < 1), and by Eq. (A.12) 
after saturation, with the temperature T and pressure p continuous at saturation (RH = 1). In the 
present approximation, none of the heat is shared with any condensate. Furthermore, even if the 
condensate is carried along with the gas and exerts no interphase drag on the gas, the condensate is 
regarded as unavailable for evaporative cooling, so reversal of the process is excluded. 

A.3 Thermohydrostatics for a Model of a Vortex in Moist Air 

At the periphery of the vortex, i.e., for r = r e (» r Q ), profiles for T(p), and for Td(p)— or 
RH(p), are taken as given, often in graphical form (e.g.. Fig. 1), typically from a sounding of the 
atmosphere in “close” proximity, temporally and spatially, to a severe tornado. From 
straightforward use of Eqs. (A.1)-(A.8), one may obtain all the dependent variables p, pv, p a , Pv> 
pa, p, T, and H as a function of altitude z, where z = 0 is approximately ground level. In practice, 
we pursue the dependent variables over a range slightly in excess of 0 < z < zj, where H (zj) = 
H(0); for the convectively unstable sounding of interest, in which a pronounced midtropospheric 
minimum in H(z) exists (see Fig. 13), zj = 0(10 km). Local anomalies may arise, such that H(z*) 
= H(0) for 0 < z* « 10 km; such anomalies are ignored in that the sounding is still regarded as 
convectively unstable overall (even if locally stable), and the value H(z*) is just an “ordinary” point 
of the profile. We shall henceforth affix subscript e to indicate that a thermodynamic-variable 
profile pertains to the vortex periphery. 

For the moist-adiabatic locus of thermodynamic states holding for region HI, Eq. (A. 10b) is 
used, with the constants (Y v ) re f, T re f and p re f associated with (Y v )eo, Teo, and peo, respectively, 
where the latter trio are ground - level values holding at the periphery : Eqs. (A.1)-(A.7) are used to 
associate p a , p v , p, Pa, and p v with p, T, and (Y v ) re f, for RH < 1. Upon satisfaction of the 
condition RH = 1, we use Eq. (A.12) to find T for all higher values of p; in this regime, p v = P(T), 
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and the other thermodynamic variables follow as before. Eventually a pair (T,p) on the moist 
adiabat matches a pair (T e ,p e ) of the sounding; the height at the periphery associated with that 
particular matching pair (T e ,p e ) is henceforth denoted zj(~ zj), and is adopted as the height of the 
lid of the storm at all radial positions 0 < r < r e . Use of hydrostatics [Eq. (A.8)], the tabulated 
locus of moist-adiabatic states [such that p(p) is known], and the known height of the “lid” state 
permit an altitude z to be associated with each state of the moist-adiabatic locus, until z = 0 is 
reached. The ground-level moist-adiabatic-locus values p(0), T(0), p(0), etc., are henceforth 
denoted p 0 , T 0 , p 0 , etc. In general, p 0 < Peo. with [(peo - PoVPeo] = 0(0.01). 


Use of hydrostatics, and of Eq. (A. 10b) with (Y v )ref = 0, T re f = T e (zx), and pref - Pe( z T)» 
gives the dry-adiabatic locus of states holding in the eye (region IV); 


7 = 7t 1 LzJ p( z t) 

MY- l)p( z T) 


p eye (z)l(T-lh 

. p( z t) . 


(A. 13) 


where the subscript gys is introduced to distinguish thermodynamic variables holding in a fully 
inserted, moisture-free eye. At z = 0, Eq. (A. 13) typically yields a value for Peye(O) such that 
[(Peo " Peye(0))/peo] = 0(0.1). 

An implication of the above treatment, in terms of Fig. 3, is that higher-altitude fluid elements 
in the potential vortex, although initially characterized by H < Heo, during downdrift into the 
surface layer attain the state H = Heo. the ground-level value of the total static enthalpy at the 
periphery, owing to enthalpy transfer across the ground plane; the total static enthalpy of these 
fluid elements remains at the value Heo during inflow in the surface layer and ascent in the eyewall. 
The value H = Heo also characterizes all fluid elements in the eye. If no such cross-boundary 
enthalpy transfer occurs, the value of Heo decreases slowly in time. 

A.4. The Thermohydrostatically Estimated Peak Swirl Speed 

In the potential vortex, by conservation of radial momentum. 


3p 

dr 


P 


Y outer 

r 


(A. 14) 


We apply Eq. (A. 14) near the bottom of the vortex depicted in Fig. 12, and adopt p * p 0 - For a 
one-cell vortex (i.e., in the absence of region IV), we adopt the modified Rankine combined-vortex 
model: 


V 0 utei( r ) = 


Vo(r/r 0 ) , 
VcyMr)", 


0 < r < r 0 
r 0 < r < r e , 


(A. 15) 
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(A. 16) 


with r e — > oo. Substitution of Eq. (A. 15) in Eq. (A. 14) gives, upon integration. 


V 0 = 


2n ) 

Peo “ Po 

jn+ l) 

po J 


1/2 


For a two-cell vortex (i.e., in the presence of region IV), we adopt: 


Voutei(r) — 


l°- 

\Vd(ro/r)" , 


0 < r < r c 
r 0 ^ r < «> , 


(A. 17) 


so 


Vo 


(2n) 


Peo ' Peye(0) 


1/2 


(A. 18) 


In general, in Eqs. (A. 16) and (A. 18), we limit attention to the case n — 1. 

For the sounding of Fig. 1, by use of the formulae presented in this appendix, we compute 
results presented in Figs. 2 and 13-15. We note that at z = zr = 1 1-6 km, the pressure is about 
21.9 kPa; the temperature, about 220 K; the density, about 0.348 kg/m 3 ; and the total static 
temperature, H/c p , about 333 K. The pressure at the base of the sounding, peo, is about 10 5 Pa; 
the pressure at the base of the eyewall, p 0 , is about 98.7 kPa; the pressure at the base of the eye, 
peye(0), is about 93.6 kPa. Since the density at the base of the eyewall, p 0 , is 1.12 kg/m 3 , we 
compute, for n = 1, that V G - 33.4 m/s for a one-cell vortex and V D = 106.9 m/s for a two-cell 
vortex. Incidentally, saturation in the eyewall occurs at z = 2.15 km, where p = 77.2 kPa, p v ~ 
1.39 kPa, T = 285 K, and p « 0.937 kg/m 3 . 
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APPENDIX B. ANALYSIS OF THE NEAR-GROUND WIND HELD FOR AN INTENSE 

POTENTIAL VORTEX 


B.l Introduction 

The velocity field in the near-ground inflow layer (for an incompressible constant-molecular- 
viscosity fluid) is examined over a closed domain of finite radial extent r e , where (say) r e = 

10 r 0 (= Ro), r 0 being the (given) radius of boundary-layer separation. While the basic results and 
flowfield structure for the surface boundary layer are long established (Burggraf et al. 1971; 

Carrier 1971; McWilliams 1971), interest persists (e.g., Phillips and Khoo 1987). Here, 
convenient approximate expressions, based upon appropriately modified “local similarity” 
solutions, are presented for the near-ground flow under the high-speed portion of a potential 
vortex. Conventional eddy-viscosity models for the swirling boundary layer also have been treated 
in detail (Carrier and Fendell 1977). 

In the high-speed portion of the potential vortex, the inward radial pressure gradient is 
virtually in equilibrium with the centrifugal force (cyclostrophic balance); there is virtually no radial 
inflow in the well -developed vortex, and there is only a weak downdraft into the surface inflow 
layer. In the boundary layer, the no-slip boundary conditions retards the centrifugal force, such 
that a “favorable” pressure gradient drives an appreciable radial influx of swirling air toward the 
axis of symmetry. Surface friction dissipates some of the angular momentum; in fact, in the 
boundary layer, well described by the parabolic approximations to the equations for the 
conservation of radial and azimuthal momentum, the radial velocity component is comparable to, 
and locally even exceeds, the swirl. Once the radial and azimuthal velocity components are 
known, the (relatively small) axial velocity component is obtained from continuity. 

The overall boundary layer grows monotonically, but rather modestly, in thickness with 
decreasing radial distance from the axis of symmetry. At larger radial distances from the axis of 
symmetry, friction plays a role across the entire vertical extent of the near-surface inflow layer. 
However, closer to the axis, the bulk of the radial inflow is described by inviscid analysis to good 
approximation; only in a thin sublayer immediately adjacent to the wall does friction alter the radial 
velocity component significantly, from a near-wall maximum value to zero at the wall (for 
enforcement of the no-slip boundary condition). The radial velocity component has a larger 
maximum at smaller radial positions, a maximum that can become as large as the swirl speed above 
the surface inflow layer at the same radial distance from the axis. On the other hand, the azimuthal- 
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velocity-component profile is shaped significantly by friction across the entire extent of the surface 
inflow layer, at all radial positions. At any fixed altitude within the inflow layer, the circulation 
becomes smaller at smaller radial positions. At a fixed radial distance from the axis, the swirl 
velocity component decreases monotonically with decreasing altitude above the ground, from its 
potential-vortex value to zero at the outer edge of the frictional sublayer. Nonoscillatory vertical 
profiles hold for the radial and axial velocity components, and assertions (e.g., Lewellen 1976) 
that the swirl within the inflow layer may appreciably exceed its asymptotic value outside the 
inflow layer, at a fixed radial position, are regarded as unlikely. 

The thinning of the frictional sublayer with decreasing radial distance from the axis is 
consistent with the generally expected behavior of a shear layer that is subject to a favorable (i.e., 
accelerating) radial pressure gradient; the frictional sublayer is subject to a more rapid radial 
velocity at its outer edge for smaller radial distances. 

In what follows we use these statements to obtain results via numerical integration of very 
simple ordinary differential equations, rather than via numerical integration of the partial- 
differential equations (Prahlad and Head 1976; Chi 1977; Shakespeare and Levy 1982). 

B.2 Formulation 

Under previously discussed approximation, if in cylindrical polar coordinates (r, 0, z) the 
corresponding velocity components are denoted (u, v, w), the conservation of mass, radial 
momentum, and angular momentum may be expressed as follows (with subscripts r, z denoting 
partial differentiation): 


(ru) r + (rw) z = 0 ; (B.l) 

(rVp.(rvt--(r,)» .[ 

T 2 

[v(r v)J z - u(r v)t - w(r v^ = 0 ; 

the boundary conditions are 

z— >«>; v— >r/r;u— >0; 

z = 0: u = v = w = 0; 

r = r e : u = 0, v specified (= I7r c ) . 

Here, r e is an inferred const., and signifies the radial extent (“edge”) of the flow domain. It is 
anticipated that the formulation is to be applied only in the range r Q < r < r e , where r 0 is the (given, 


(B.3) 

(B.4) 

(B.5) 

(B.6) 
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finite) radius of boundary-layer separation. Also, in the potential vortex, T is effectively given, 
where (with v ou ter = V, for brevity) 


T = r V = r G V 0 , const , (B.7) 

and V 0 is inferred from the given peripheral sounding (Appendix A). It may be noted that as a 
consequence of standard boundary-layer approximations, the cyclostrophic balance has been used 
to identify the radial pressure gradient in Eq. (B.2): 


1 

P 



(B.8) 


B.3 Constant-Viscosity Solution 

If the kinematic viscosity v is const., then solution to Eqs. (B.l) - (B.7) is sought in the form 
[the function P(r) is to be distinguished from the parameter P, introduced in Eq. (12)] 


r u = -r P(r) g'(q) , r v = Tf(ti) , w = W(t|) , (B .9) 

where 


n = 




(B.10) 


The functions s(r) and b(r) are to be identified in solution. One does not expect Eqs. (B.9) and 
(B.10) to give an exact solution, but to capture the essence of the functional behavior in a very 
useful form. 


Substitution of Eqs. (B.9) - (B.10) in Eq. (B.l) gives 

Ps' 


w^(vrr[^p- g '-^ng- 


= (v r) : 
=(*r) 


, 1/2 


:l/2 


^P'g'- 


Ps' 

2 r s 1/2 


(o g' - g^ 


1/2 


S 1 / 2 P' ps' ps' 

— g * -r^-rjz n g 
r 2 r s 1/2 ' 2 r s 1/2 


(B.l 1) 


where, by Eq. (B.5), 


g(0) = 0. (B.l 2) 

Substitution of Eq. (B.9) - (B.10) in Eqs. (B.2) and (B.3) gives, respectively. 
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(B.13) 



r 


2 






= 0. 


(B.14) 


The terms in Eq. (B.13) derive, from left to right, from similarly positioned terms in Eq. (B.2). 
For self-consistency, one must take 


*P' 



(B.15) 


say; also, it is necessary, but not sufficient, to take 



e ' 2 = i - f 2 . 

(B.16) 


PP'_ i-P 2 

r r 2 

(B.17) 

From Eq. (B.6), p(r e ) = 0; hence, from Eq. (B.17), 



p 2 = r * ‘ r2 . 

ri 

(B.18) 

Equation Eq. (B.15) may be written 



(p 2 s)' 

= - 2 r p => s = | r e (r 2 - r 2 ) 1/2 

(B.19) 


by use of Eq. (B.18). 

The term in square brackets in Eq. (B.13) is 0(Pr2/s), while the other terms are 0(F^/r-); 
substitution of Eqs. (B.18) and (B. 19) shows that the ratio of the bracketed to unbracketed terms is 

o[(3r2)/(2r?)], so that, for r « r e , the bracketed term is negligible. The neglect of the bracketed 

term gives an inaccurate result where the motion is slow [i.e., where v = O(v max /10), 
u < O(v max /10)], and an accurate result where the motion is fast. From Eq. (B.19), s 1/2 ~ 

(r e - r) 1/4 , as r -» r e ; this is consistent with the more meticulous analysis of Stewartson, 

Burggraf, and Belcher (1971). The flow structure remains less evident in other approaches (e.g., 
Kuo 1971, 1982). 

The radial momentum equation for r| » r 2 , so that p -> 1, becomes a statement that the left- 
hand side of Eq. (B.2) vanishes; i.e., V 2 = u 2 + v 2 , such that the radial velocity component is 
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inviscidly controlled. The two-point nonlinear third-order boundary- value problem, from Eqs. 
(B. 13) - (B.16), becomes 


f' + gf = 0; 

(B.20) 

g' = (l -f 2 ) 1 ' 2 ; 

(B.21) 

f(o) = g(o) = 0, f (°°) -> 1 . 

(B.22) 


The boundary conditions follow from Eqs. (B.4), (B.5), (B.7), (B.9), (B.10), and (B.12). The 
problem is readily solved numerically by shooting; one finds 

f(0) = 0.7456, g(~) = 1.941 . (B.23) 

Results are graphed in Figs. 17 and 18, where it is evident that f monotonically increases as T| 
increases (so g' monotonically decreases). Note that g'(0) = 1 (cf. Rotunno 1980). 

Of course, the reduction to third order results from dropping of the frictional term in the 
conservation of radial momentum Eq. (B.3), which has been taken in the inviscid form 
Eq. (B.16). Enforcement of all of Eq. (B.5), specifically g'(0) = 0 so that u(x, z = 0) = 0, 
requires restoration of the frictional term in a thin near-wall sublayer of the inflow layer. 

For completeness, a brief sketch of the near-wall sublayer, in which u falls to zero, is 
developed now. If [with the dimensionless function h(o) introduced here to be distinguished from 

the eyewall thickness h(z), introduced in the main text] 

r u = . r Mr) h'(o) ,<J=^(r) 1,2 .|5=l- (r/r 'f • < B24) 

then from continuity 

(r w)o = m|^j 1/2 [p* T h' - |J T h"(— ) o] , (B.25) 


or 


r w) = m^) 1/2 P' T h - p lfe)(« h' - hj 


(B.26) 


From the radial momentum equation. 


AW 

lm 2 ^ 


^V + 


p(p 


^^] hh - = x.(g h - 2 + (i)^. 


(B.27) 
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where the last term involving f 2 is dropped as negligible. The function m is assigned by 
demanding 

ft 2 m ' _ ! 

rm r 2 (B.28) 

If attention is concentrated on (r/r e ) « 1, such that p — » 1, P' —> o(l), then 

m = r , (B.29) 

and 

- h'" + (h h^ = 1 , (B.30) 

where all neglected terms are O (r 2 /r|) or smaller. The boundary conditions are 

h(0) = 0 , h'(0) = 0 , hioo) -> 1 ; (B.31) 

these enforce w(x,z = 0) = 0, u(x,z = 0) = 0, and matching of the sublayer and outer-layer 
solutions, respectively. One arrives at a Riccati equation 

-h / + ^- = ^--h"(0)a,h(o->oo)->a, (B.32) 

where the boundary conditions conditions have been used to assign a constant of integration. The 
results of Burggraf, Stewartson, and Belcher (1971) assure that a solution exists, and for present 
purposes 

h(a) = o + exp{-o) - 1 => h'(a) = 1 - exp (-a) , (B .33) 

captures its behavior adequately (see Fig. 19). Hence, a uniformly valid expression for r 2 « r| is 

r u = - T |3{r) [h '(a) + h'(r») - 1] , (B .34) 

where h'(cr) is given by Eq. (B.32), or , facilely, by Eq. (B.33), and h'(r|) is given by Eqs. 

(B.20) - (B.22). 

This modification to encompass the sublayer is regarded here as a detail from an engineering 
point of view; the solution to Eqs. (B.20) - (B.22) would normally furnish an adequate 
approximation by itself. However, completion of the surface-inflow-layer analysis does 
emphasize the compatibility of separated-inflow-layer properties with the properties of contiguous 
eye air on one side of the annulus, and with the properties of the contiguous potential-vortex fluid 
on the other side. The fluid in the nonrotating sublayer of the inflow becomes contiguous to the 


23 


R5 92.019 



nonrotating, slowly recirculating air of the eye. The rapidly rotating fluid of the outer layer of the 
inflow becomes contiguous to the rapidly rotating air of the potential vortex, and exerts very little 
streamwise shear on the potential- vortex air. 

B.4 Pertinence of the Surface-Inflow-Layer Solution 

The quantities q 0 and v 0 in Eqs. (4) and (8) are transversely averag ed values (of the 
streamwise and swirl velocity components, respectively) holding at z = 0, the bottom of the 
eyewall, where the height z = 0 is after the turnaround. If the kinetic energy is conserved through 
the turnaround, then the relation V 2 = u 2 + v 2 , given just above Eq. (B.20), suggests the relation 
Vg = q 2 + vg, where, at separation of the surface-inflow layer, the streamwise-velocity component 
q2 is identical with the radial inflow u^. We adopt, at a fixed value of r « r e , a flux- weighted 
transverse average of the radial-inflow velocity component u and of the swirl component v, for the 
purpose of establishing the value of the ratio (q 0 /v 0 ): 


9 ° =. 

Vn 


2rtrp j u 2 dz j 


(l-f^dii 


27crp j |u| v dz j f(l-f 2 ) ,/2 dq 

Jo Jo 


(B.35) 


by use of Eqs. (B.9), (B.10), and (B.21). Attention is confined entirely to the major, inviscid 
portion of the inflow layer, described implicitly by the solution of Eqs. (B.20) - (B.22). 

Substitution of the numerically derived solution for f(T|) in Eq. (B.35) gives (qo/v 0 ) « 1.07, or, 
approximately, q 0 = v 0 . If this relation between q 0 and v 0 continues to hold through the 
turnaround, Eq. (8) may be written as V£ = Qo + vg = 2qg = 2v 2 More generally, the streamwise 
velocity component may be decremented owing to a vortex breakdown in the turnaround, whereas • 
the conservation of angular momentum plausibly holds across the turnaround, so v 0 remains 
unaltered. Thus, in general, 



where M(> 2 1/2 ) is a factor to be assigned for parametric investigation in view of the limited 
understanding of the “breakdown” of an annular column with swirling flow. 

A motivation for regarding an abrupt deceleration of the streamwise flow in the swirling- 
updraft annulus as plausible is provided by the dust striations recorded in certain tornado 
photographs [e.g., of the June 5, 1966 tornado at Enid, OK (see Gwynne 1982, p. 55; Lugt 1983, 
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p. 192)], and by examining the results of photogammetry carried out for a few tornadoes (Grazulis 
1991, pp. 100- 101). Under the approximation that the interphase velocity slip between the 
observed particles and the colocated air is negligible, we note that, for the low-altitude, near-core 
flow in the tornado, the angle of the velocity vector for the fluid with respect to the ground plane 
typically tends to be appreciably less than 45°. Thus, the inequality v > q is inferred in the post- 
turnaround flow from photographic data, although v *» q at separation of the surface inflow layer, 
from the above analysis. This would be consistent with a nonentraining vortex transition (with 
conversion of mechanical energy to heat) in the turnaround region. 
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APPENDIX C. THE SELFCONSISTENT VORTICAL FIELD OUTSIDE THE EYEWALL 


Derivation of Eq. (7) involves examination of the isobars in region I, the outer potential 
vortex, holding for R(z) < r < <» (for r e — > «>). Under the cyclostrophic and hydrostatic 
approximations, appropriate for region I, 


Pr = pVouteA . (C.l) 

Pz = -Pg • (C.2) 

where the notation is standard (except that we emphasize by subscript that the relevant swirl 
pertains to the potential- vortex region only). On an isobar, 

dp = p r dr + p z dz = 0 , (C.3) 


or 


dr _ Pz _ §£ 


dz 


Pr 


¥ outer 


(C.4) 


Under the family of spin-up states specified by Eqs. (11) and (12), 


V 2 oR 2 o 



For the linear profile for f(z) given in Eqs. (13), 

VoRo -gfoj a-b zfrp)] ^ 
2r 2 b a - b h(p) 


(C.5) 


(C.6) 


where a = 1, b = (1 - f5)/zj, and 0 < p < 1, with zj obtained from Appendix A and with P 
specified as a parametric input The pressure p takes on the value for which the integration is 
performed; z(°°;p) is the height of the isobar for pressure p at the periphery of the storm; and z(r;p) 
is the height of the same isobar at some finite cylindrical-radial distance from the axis. We 
henceforth denote z(«>;p) as h(p), the altitude at which the unspun periphery has the pressure p 
according to results of Appendix A. (The tilde is present because symbol h has been preempted for 
another quantity.) If we let r = R(z), the locus of the eyewall displacement, then 


VqR^ -g 1J 1 -MM f 
2R 2 (z;p) b 1 - b h(p) J 


(C.7) 
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where the pressure p at r = R is denoted pamb- By rearrangement, we obtain Eq. (7), explicitly, 



1 -(1 - bz)ex 



(C.8) 


Upon fixing the parameter b, each pair of values for z and R implies a value for h(pamb)> then, 
by interpolation in the thermohydrostatic tables, pamb- F° r example, at z * 0, where R = Bo. 
h = ( 1 - exp[- b V^2g)])/b, so h > 0. As expected, the isobars dip groundward as the radial 
distance from the axis of rotation decreases; only for R 2 » (b Ro)/2g does h(pamb) z, 
height of the isobar at the eyewall/potential-vortex interface. Since thermohydrostatic results 
typically are not tabulated for altitudes above zj, tabulated values for h may be exhausted before 
the top of the eyewall is attained in the integration described in Section 3. 

If, instead of the linear relation for f(z) given by Eq. (13), we adopt the two-part expression 
(say, 0 < P 2 - 0-3) 


f(z) = 



0 < z < zj 
Zi < Z < ZT , 


(C.9) 


where (for the sake of an explicit choice) z\ may be taken to be the altitude at which the eyewall 
becomes saturated. Substitution of Eq. (C.9) in Eq. (C.5) gives, for r — R(z), 

Rpamb)=z + ^^,0<z<h£z i; (C.lOa) 

2 g R 2 


h(Pamb) ~ Z J _ 1 \ ~ Q 

Zj* Z\ 


pi/2 1+Q ’ 


where. 


if 5 2 = p2 /2 [(z-ziX z T-zi)] • 


Q = 


exp 

1 - S 2 

,1+52 


2pi' 2 

'V^R 2 

Zt-Zj 

_zgR 2 


(z-Zl) 


exp 


Pi' 2 V 2 Ri 


g(zj - Zi) R 2 J 


Z < Zj < h < Zj 
Zj < Z < h < ZJ 


(C.lOb) 


(C.lOc) 
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APPENDIX D. THE TURNAROUND 


We consider the very-low-level portion of region in, in which the surface-inflow layer 
separates and “overshoots” its equilibrium position, owing to dynamic forces; consequently, in the 
“overshoot”, under conservation of angular momentum, the local swirl speed may modestly exceed 
the thermohydrostatically computed peak value Vq. Eventually, restorative forces result in an 
“overcompensatory” locus for the lower eyewall, such that the local swirl speed may “undershoot” 
its equilibrium value. In fact, an endless cycle of overshooting and undershooting behavior is 
formally obtained as the solution for the formulation (for the turnaround), which is translationally 
invariant in the independent variable (the altitude z). The limit-cycle behavior is of purely academic 
interest because: (1) the translationally invariant formulation is inadequate (because, inter alia, 
change of the eye pressure and of the gravitational potential energy with altitude cannot be 
justifiably ignored), and, far more importantly, (2) the flow is likely to be unstable to breakdown 
in view of the centrifugal forces acting in a flow with curved streamlines. Vortex-breakdown 
theory (Hall 1972; Leibovich 1976) pertains to axis-enveloping swirling flow (with a stagnation 
point on the axis), often (but not always) confined in a tube; the annular geometry of relevance here 
is largely uninvestigated. Nevertheless, although we can provide no quantification, the existence 
of a breakdown is highly plausible. We proceed with an analysis of the turnaround to confirm that 
the equations and boundary conditions yield solutions consistent with the just-described 
anticipations. 

Equations (1) - (3) and (6) are unaltered from those given earlier, but Eq. (7) is modified 
because we take p = p 0 in the very low portion of the potential vortex, and we adopt the 
entrainment-free (ct E = 0) form of Eq. (1). Then, if Vouter denotes the swirl speed in the potential 

vortex. 



outer 

r 


= Po i ^ 


(D.l) 


hence, in the potential vortex, 


p = _PoV|jRlj + const. (0.2) 

At the separating-layer/potential-vortex interface, r -» R and p — » Pamb> by definition: 


Pamb = -^2^(|f) + const. 


(D.3) 
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If the frictional layer separates smoothly (the particular case to be examined here), since the 
pressure is invariant, to lowest order of approximation, across a thin near-wall layer under the 
classical arguments of Prandtl, then Pamb = Peye for R = Ro (see Fig. 3). Under this assignment 
for the constant of integration, 


_ p 0 Vo 

Pamb * Peye — j 




(D.4) 


where the last equality arises from the fact that Eq. (9) may be rewritten (Appendix B) as 

V£ = q§ + vg - 2ql = 2 v§ , 


(D.5) 


prior to the occurrence of any possible breakdown phenomenon. The above derivation implies that 
the potential-vortex behavior describes the swirl across the separating inflow layer. 

Bernoulli’s equation [Eq. (4)] becomes, upon neglect of the gravitational potential energy, 
under Eqs. (6), (D.4) and (D.5), 


v 2 + q 2 I dp' _ q§ + v| 

2 + l p( P 'f 2 

J Peye 


where p = Peye + 


PoqS 

2 



(D.6) 


from Eqs. (6) and (D.4). Since v 2 = qg Ri /R 2 , and since p(p) is tabulated from preliminary 
hydrostatics, Eq. (D.6) provides an expression for (q/q 0 ) as a function of (R /Ro)- 

We adopt the following nondimensionalization and parameter definition: 


R 

q 




a s Ro/h 0 , T = 


2 Peye 

Po 9o 


(D.7) 


so 


p h = (q R) 1 , v = (r)" 1 , and 

± „r(p. i)[i + (RiF]‘' 2 4R;T]j , 

r 2 ^R;t] ' r 


(D.8) 
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where 



(D.9) 


/■1+X 


FM» 


dp - 

p(?) 


(D.10) 


Since we anticipate that Rj — » ±°° within the domain of integration, it is locally convenient to 
interchange the roles of dependent and independent variables, so R(z) -» ^R): 


R z - (zr ) 1 » Rzz 


ZrR 

(Zr ) 2 


(D.ll) 


Hence, 

m = JjtiSLj.l ± orIr 2 - i)[i ♦(^''^ROr]) 1 ' 2 • <d.i2) 

RR r 2 ^r ; t]I R 1 

The symbol ± before the term containing a square root implies that we may choose the physically 
pertinent branch of the root, i.e., the physically pertinent curvature. 

Starting conditions, holding at the position of the surface-inflow-layer separation, are taken to 
be those for a smooth separation: 

z = 0 ,R= 1 ,Zr = 0, (D.13) 

such that Eq. (D. 12) is convenient for the initial incursion (R < l) of the separated layer, until Zr 
becomes large. Then, with continuity of the values of z, R, and Zr[= (Rz)' 1 ] at a convement point 
of transition (say, |z^ = 1), we adopt Eq. (D.8) to proceed to R = R min and continue on until 
R — > 1 and R^ — > °°, so we revert to the use of Eq. (D.12) to begin the subsequent excursion 
(R > l), etc. If we denote Eq. (D.8) as formulation I and Eq. (D.12), formulation II, then Fig. 20 
indicates which formulation is convenient for various portions of the incursion/excursion cycle; 

Fig. 20 also indicates which trial-and-error choice of sign for the square root is physically 
appropriate, where the choice of sign in formulation I may be switched across a point at which 
|zp>! -4 co, i.e., across R = R min « 1) and R = R ma x (> D, and, correspondingly, Ae choice of sign 
in formulation II may be switched across a point at which |R^ -* i.e., across R = 1. 
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Results are presented in Fig. 21 for the physically interesting parameter values ot 5 and 10, 

T = 29.4, under the following curve fit to thermohydrostatic tabulations for the lower portions of 
the moist adiabat for a specific sounding (Figs. 1 and 2) of a tornado-proximity atmosphere 
(Appendix A): 

p(p) * 0.0654 + 1.066 p - 0.1663 p 2 (D.14) 

over the range of p of interest, explicitly, p = 0(1). Since v = (S)" 1 , in view of Eqs. (D.5) and 
(D.7), the value of (2 1/2 Rmin) -1 gives the ratio of the peak swirl to the nominal, 
thermohydrostatically computed peak swirl, Vq, for a specific (but typical) case. This very local 
increase in the swirl over the value V D is termed the “overshoot”. All the dependent variables may 
be straightforwardly deduced from knowledge of R(z), but these details are not pursued. We 
anticipate that a breakdown occurs on the first occasion in the turnaround for which an incursion of 
the separated inflow layer undergoes transition to an excursion, but we are highly uncertain about 

the details. 
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RIM. 92. 0266. 10 


Fig. 1. Convectively unstable sounding (with a moderately large midtropot (P h f n c 
total static enthalpy), measured in the vicinity of a tornado at Jackson, MS on Ap 
2300Z. The static temperature is T; the dew-point temperature, Tj; the pressure, p. 
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ground level (diamonds). 
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RIM 92.0266.07 


Fig. 3. Schematic of a postulated four-part model of the structure of a quasisteady two-cell vortex, 
of axial extent zj (the “lid” of the storm). Region I is the outer potential vortex; II, the near-ground 
inflow layer; HI, the “eyewall”, the lowest portion of which is termed the turnaround; and IV, the 

“eye”. The dotted curve is a sketched isobar, which is at altitude h at the periphery, r -4 r e , and at 

lower altitude at the “eyewall’Vpotential-vortex interface, r = R(z), with R(0) = Rq and r e » Ro- 
The arrows schematically indicate the magnitude and direction of the secondary (i.e., radial and 
axial) flow. 


37 


R5 92.019 



h (m) 



Fig. 4. For the nominal case [defined by the sounding of Fig. 1 and the parameter assignments of 
Eq. (17)], the radial displacement of the eyewall, R(z), and the thickness of the eyewall, h(z). 
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h(km) 



Pamb 

RIM. 92 0266 21 


Fig. 6. For the nominal case, the pressure at the potential-vortex/eyewall interface, PambUX and 
the height at the periphery at which the pressure Pamb occurs, h(z). 
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T (s) xIO ' 3 



Fig. 7. For the nominal case, the mass in the eye as a function of altitude, m eye (z), and the filling 
time T(z)= m e ye(z)/m e yewaii> where m ey ewaii is defined by Eq. (14). 
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R (km) 

R1M.92.0266.23 


Fig. 8. The radial displacement of the eyewall, R(z), for several values of the parameter (3, related 
to the vertical distribution of the angular momentum at the periphery, where P = 1 is an axially 
invariant distribution and p = 0 is strongly stratified. All other parameters and functions are held at 
their nominal values (and p = 0.99 is nominal). 
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Fig. 9. The eyewall thickness, h(z), to complement the results of Fig. 8. 
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Fig. 10. The radial displacement of the eyewall R(z), for three vdues of the initial ^^ ess ^ 
all other narameters and functions being held at then- nominal values (with ho - 35 m being 
nominal ) P Of course, the initial curvature e [Eq. (18)] varies with h G , since the ratio (Ro/ho) is 

fixed at its nominal value, 10. 
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z 



RIM. 92. 02 66. 04 


Fie. 12. A simplified schematic of the structure of an axisymmetric two-cell vortex convenient for 
thermohydrostatic computation. The surface-friction layer II and the turnaround (Fig. 3) have 
been omitted; in the above sketch, the surface frictional layer does not separate smoothly, and the 
pressure at the base of the eyewall (region HI) appreciably exceeds that at the base of the eye 
(region IV). For a less intense, one-cell vortex, no eye is present, and r 0 -» 0, so region 111 (no 
longer termed an eyewall) envelopes the axis of symmetry. 
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Fig. 13. The total-static-enthalpy profile H(z) associated with three columns of fluid for the model 
of Fig. 12, with the sounding of Fig. 1 holding at the periphery. 
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Fig. 14. The pressure profile p(z) associated with three columns of fluid for the model of Fig. 12, 
with the sounding of Fig. 1 holding at the periphery. 
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= ig 15 The temperature profile T(z) associated with tttree columns of noid for the model of 
FI 12, with the sounding of Fig. 1 holding at the periphery. 
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R1M.92.0266.16 


Fig. 16. The density profile p(z) associated with three columns of fluid for the model of Fig. 12, 
with the sounding of Fig. 1 holding at the periphery. 
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Fie 17 For spatially constant viscosity, results related to the radial and axiaJ , ve J^^' ^f™n° nentS 
for the outer preponderant portion of the boundary layer under the lugh-speed portion of an 
Messed vS The similarly independent variable n is large for large distances normal to the 
wi£l or for small radial distances from the outer edge of the vortex system. 
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Fig. 18. Asa complement to Fig. 17, results related to the swirl-velocity component. 
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Fig. 19- For spatially iXw" to *eihm, 

and approximate analytic (solid) results for the run f impressed swirl. The similarity 

SSSM3S "Xe dis^c'i P norm P al to toe wall or for small radial distances 

from the axis of rotation. 
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pig 20 For numerical integration of the turnaround, a convenient (but nones^ntial) eight- 
segment SdivS is considered, with a segment end-point introduced wherever 

to = 0, Rz = 0, or |Rzl = 1. Equation I is Eq. (D.8); Eq^Il is Eq (D.12), enva 
bj an interchange of the roles of the independent and dep^dent v|iable Si Th^ s | n 
indicates the sign locally adopted for the last turn m Eq. (D.8) or Eq. (D.12), as appropnaie 
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Fig. 21. The dimensionless radial displacement, R®, of a smoothly sep^tM^rfaa-mnow 
layer over one cycle (incursion, then excursion, with increasing attitude z) of the turnaround. As 
the geometric ratio a increases, the minimum value of R increases, so the local peak swirl exceeds 
the thermohydrostatically estimated peak value by a smaller amount 
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